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We explore ballistic regime quantum transport characteristics of oxide-embedded crossing 
and kinked silicon nanowires (NWs) within a large-scale empirical pseudopotential electronic 
structure framework, coupled to the Kubo-Greenwood transport analysis. A real-space 
wave function study is undertaken and the outcomes are interpreted together with the 
findings of ballistic transport calculations. This reveals that ballistic transport edge lies 
tens to hundreds of millielectron volts above the lowest unoccupied molecular orbital, with 
a substantial number of localized states appearing in between, as well as above the former. 

We show that these localized states are not due to the oxide interface, but rather core 
silicon-derived. They manifest the wave nature of electrons brought to foreground by the 
reflections originating from NW junctions and bends. Hence, we show that the crossings 
and kinks of even ultraclean Si NWs possess a conduction band tail without a recourse to 
atomistic disorder. 

PACS numbers: 73.21.Hb, 73.63.Nm, 72.15.Rn 

I. INTRODUCTION 

As the fabrication techniques for silicon nanowire (NW) 
polymorphs are relentlessly advancing, 1 8 the latest ac¬ 
complishment has been their assembly into higher dimen¬ 
sional complex architectures. With this aspiration, a new 
class of nanostructures emerged, including L-shaped 9,10 
and kinked 9 as well as hyper-branching 11,12 and networking 
Si NWs. 13 Among these, kinked NWs (KNWs) in partic¬ 
ular have received special interest owing to their rich ge¬ 
ometries including U-, V-, and W-shaped varieties. 10,14 ” 16 
They have potential application areas ranging from bio¬ 
logical and chemical sensing, nanoscale photon detection 
to three-dimensional recording from within living cells and 
tissue. 17-19 

These NW-based superstructures in their ultraclean 
forms can be regarded as the electron waveguide com¬ 
ponents making up a quantum network, once theoreti¬ 
cally envisioned. 20 As a matter of fact, the field has a 
rich background. 21 The dawn of the ballistic transport 
era in low-dimensional semiconductors can be marked by 
the quantum interference transistor proposal of Datta and 
coworkers that was based on exerting a phase difference 
between the two electron channels. 22 This was followed 
by experimental reports on GaAs/AlGaAs heterojunctions 
demonstrating resistive increase when the current path en¬ 
counters a bend, 23 and also resistance oscillations on T- 
shaped electron channel when the stub length is electri¬ 
cally tuned. 24 More recently, as observed in a 500 nm-long 
silicon NW, quantum interference effects cause an oscilla¬ 
tory pattern in conductance as well as a shift in conduction 
threshold depending on the quantization energy within the 
wires. 25 Despite the intense efforts, due to imposing chal¬ 
lenges in avoiding diffusive transport channels, the elicita- 
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tion of experimental fingerprints endorsing quantum inter¬ 
ference effects remain still very subtle. 

Without any doubt, any theoretical insight would be in¬ 
valuable to better understand the rich physics and appli¬ 
cation prospects of these structures such as silicon-based 
crossing NWs (CNWs) and KNWs. If one sets the objec¬ 
tive to characterize the full mosaic of localized to delocal¬ 
ized states, then it becomes imperative to acquire atom¬ 
istic resolution, and therefore the methodology needs to 
surpass the effective mass, k • p, or transfer matrix-like 
continuum models. 15,26 28 On the other hand, due to large 
system sizes in CNWs and KNWs, on the order of ten thou¬ 
sands of atoms including the silica embedding host matrix, 
they are not amenable by the current state-of-the-art ab ini¬ 
tio techniques. 29 Because of these constraints, only a small 
number of realistic-size atomistic NW studies exist, 30-32 
leaving the CNWs and KNWs as an unchartered territory. 

This work aims to provide a computational exploration 
of the full scale, from localized to semi-extended up to con¬ 
ducting characteristics of electronic states in Si CNWs and 
KNWs. We employ an empirical pseudopotential-based 
atomistic electronic structure solved using the linear com¬ 
bination of bulk bands (LCBB) method. 33,34 Very recently, 
we have successfully applied this technique on the bandgaps 
and band edge alignments of Si NWs and NW networks. 35 
In conjunction with the quantum ballistic transport cal¬ 
culations utilizing the Kubo-Greenwood (KG) formula, we 
undertake an isosurface analysis of wave function distribu¬ 
tions for the CNWs and KNWs. Our main finding is the 
identification of a conduction band tail of up to several 
hundred millielectron volts in span. Unlike the well-known 
Urbach tail associated with an atomistic disorder, 36,37 it 
stems from the electron wave localization within the junc¬ 
tion and kink sections of ultraclean Si CNWs and KNWs. 

The temperature range relevant to this ballistic regime 
is a few Kelvins, even though we shall be presenting zero- 
temperature results. Another noteworthy point is the 
strain which can in general play a significant role on the 
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electronic properties of Si NWs . 38 40 Yet, a uniform strain 
(radial or axial, compressive or tensile) retains the exten¬ 
sions of states throughout the NW and does not lead to 
localized states . 38,41 In contrast, inhomogeneous strain can 
promote the localization behavior in NWs; in the same vein, 
interface roughness can induce disorder-originated localized 
states . 31,42,43 In this work, targeting atomistic disorder-free 
localization, the strain, as well as the interface relaxation 
effects, are left out of the scope, to unambiguously account 
for the effects of wire topology. Nevertheless, including the 
strain and interface disorder is expected to further intensify 
and enforce the present localization behavior. 

The paper is organized as follows: in Sec. II we present 
the general theoretical framework of our electronic struc¬ 
ture and transport analysis. In Sec. Ill we provide our re¬ 
sults on CNWs and KNWs, and conclude in Sec. IV. In the 
interest of a lucid presentation, we defer some technicali¬ 
ties to two Appendixes. The first one provides the essential 
details on the LCBB technique, followed by the details on 
the atomic pseudopotentials employed in this work. 


II. THEORY 

For the electronic structure calculations of large-scale 
nanostructure systems the LCBB method has been 
introduced , 33,34 with the distinction to solve the atom¬ 
istic pseudopotential-based Hamiltonian using a basis set 
formed by the bulk bands of the constituent materials. The 
virtue of the method lies in the nanostructure sizes acces¬ 
sible with a reasonable computational budget. Within a 
3-dimensional (3D) supercell approach the complex NW 
structures are embedded into an oxide matrix to passivate 
the dangling bonds and to prevent the interactions between 
the structures and their images in neighboring cells. The 
nanostructures considered in this work contain about 1500 
core Si atoms, and including the host medium the super¬ 
cells contain more than 10 000 atoms. Perfect crystalline 
order is considered as no structural relaxation is performed. 
Being interested in conduction states, we neglect the spin- 
orbit interaction which is mainly influential on the valence 
states. Details of the LCBB formalism, description of oxide 
passivation and the employed pseudopotential form factors 
can be found in Appendixes A and B. 

To study the linear-response regime ballistic quantum 
transport characteristics of the nanostructures, the KG- 
formula 44,45 is applied as a postprocess over their electronic 
structures. The zero-temperature dc-limit KG-formula for 
conductance is given as 

G a p{E F ) =9 s^ T ]T(n' \P a \n)(n\Pp\n') 

m o L sc n n , (1) 

x S(Ep — E n )S(Ep — E n f ), 


where cq /3 are the Cartesian indices, g s is the spin degen¬ 
eracy, m 0 is the free electron mass, and L sc is the length 
of the computational supercell along the chosen transport 
direction, P is the momentum operator and E n is the en¬ 
ergy eigenvalue of the eigenstate |n ). 46 The KG-formula 
establishes the relation between the electronic structure of 
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FIG. 1. Demonstration of the coalescence of isolated conductance 
peaks with k -grid resolution to ideal steps for a representative square 
cross-section hard-wall NW. Logarithmic scale is used for conductance 
to accommodate all the peaks. E n panels show the energy spectra for 
the associated k- grid. 


a system and its conduction properties, as a function of 
Fermi energy Ep- In our implementation, the Dirac delta 
functions of Eq. ( 1 ) are broadened by the Gaussian function 
with a typical standard deviation of 5 meV. 

It is worthwhile to discuss some technical aspects of the 
implementation of KG-formula . 47,48 Fig. 1 shows two repre¬ 
sentative cases for a free-electron model 1 nmxl nm cross- 
section hard-wall NW. On the left column, the coarse-grid 
case is displayed in which the sparse energy spectrum ends 
up with the Dirac delta function-derived isolated peaks 
both for the density of states (DOS) and conductance. The 
performance of the griding depends on AE/j, i.e., the ra¬ 
tio of the energy difference of successive k- points, AE, to 
the standard deviation of Gaussian broadening, 7 . In the 
coarse case, the ratio is about 5 at the band edge and 
further increases through the band dispersion. Since the 
energy difference between successive conductance steps is 
larger than the energy broadening, the isolated peaks of 
both DOS and conductance do not coalesce to reproduce 
the well-known analytically expected ID DOS and quan¬ 
tized conduction behavior, a signature of the ID ballistic 
transport . 21,49 Our tests indicate that as the ratio drops 
below 0 . 02 , the coalescence of isolated peaks starts to yield 
reasonable DOS and conductance shapes. On the right 
column of representative fine-grid case, by increasing the 
number of wave vector, k-points in the energy dispersion 
curve, the spectrum is made dense enough so that the ra¬ 
tio is set about 0.005 to obtain the analytical behaviors. 
Hence, this figure signifies the necessity for a dense energy 
spectrum in achieving ideal quantized conductance steps 
through the KG-formula. 

Notwithstanding, a typical LCBB calculation based on 
a 3D k-grid (see Fig. 2 (a) and also Appendix A) yields 
relatively sparse discrete energies. Obviously, the energy 
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FIG. 2. Schematic illustration of (a) 3D and (b) 2D k-grids. The 
tube indicates the alignment of considered dispersion direction. The 
planar k -grid kj_ is shifted gradually along the dispersion direction 


spectrum can be made denser (i.e., AF/y reduced) by in¬ 
creasing the number of k-points in the 3D grid (i.e., in 
the expansion basis), however it becomes computationally 
quite demanding. Instead, regarding the k-space, we con¬ 
struct a dense k^-grid over the plane perpendicular to the 
KNW axis or to the axis of one of the crossing wires of 
a CNW, and obtain the energy dispersion relation E n (k\\) 
of subbands by shifting this plane along that axis of ky 
(see Fig. 2 (b) and also Appendix A). Once a reasonable 
dispersion is obtained, by further interpolating the energy 
spectrum the KG-formula can be applied to disclose the 
conductance steps. We call this the dispersion-based ap¬ 
proach. It should be noted that in the 3D k -grid case 
which we shall name as the state-based approach, the in¬ 
terpolation is not possible since the 3D eigenvalue prob¬ 
lem yields only state energies, E n without revealing their 
/c-dependence. Thus, KG-formula inevitably produces iso¬ 
lated conductance peaks for this state-based case. 

For clean single NWs having no localized states, the 
dispersion- and state-based approaches yield the same en¬ 
ergy spectrum—once constructed with the same number 
of k-points. However, in the case of CNWs and KNWs 
the findings of these two approaches deviate: while the 
state-based approach is able to capture both the states 
having localized or extended characters, the dispersion- 
based approach finds only the extended-dispersive states. 
As we mentioned above, the state-based approach at full 
performance becomes computationally expensive, so we 
apply a two-pass approach, where we get the extended 
states through the dispersion-based approach, followed by 
a reasonable-budget state-based approach to single out the 
localized states and to mark the ballistic transport edge 
(BTE). 

III. RESULTS 
A. Crossing Nanowires 

We start our transport analysis targeting the low-lying 
conduction band edge states with the CNWs made up of 
three crossing 1.5 nm diameter Si NWs oriented along the 
(110) directions. We choose the (110) Si NWs since it is the 
major growth orientation within the sub- 10 -nm diameter 
size regime . 5 In Fig. 3, we present the isosurfaces of carrier 
densities (wave function modulus squares) of the nine low¬ 
est electron states, where the first one, CSi corresponds to 
the lowest unoccupied molecular orbital (LUMO ). 50 Among 


Three-wire Crossing Charge Densities 


| di= 1.5 nm, d2= 1.5 nm, d 3 = 1.5 nm 



FIG. 3. Isosurfaces for charge densities corresponding to the first nine 
states above the conduction band edge for the three-wire crossing with 
wire orientations along (110) directions and diameters of all 1.5 nm. 
For clarity, the surrounding oxide atoms are not shown. In this and 
following figures, LUMO energy is set to zero, and the isosurfaces 
enclose densities from 100% to 10% of their respective maxima. 


them the first six states, CS 1 -CS 6 are observed to be 
highly localized at the junction of the CNW and dimin¬ 
ishes sharply beyond the junction. Energetically, they nest 
in a window of 51meV with the individual splittings be¬ 
ing governed by the CNW as well as the atomistic point 
group symmetries as in Si nanocrystals. 51 Separated from 
this group by about 200 meV, another set of states exists, 
this time having extended wave functions spread to some 
portions of the network. As the figure demonstrates, all 
the states are Si core-derived, with no interface trapping or 
penetration to the oxide matrix. 50 

To confirm our isosurface-based interpretations, we com¬ 
pare them with those of the quantum ballistic transport 
results using the KG-formula. Specifically, we apply the 
aforementioned two-pass treatment to harness their com¬ 
plementary features to the CNW of Fig. 3. On the top 
panel of Fig. 4, the dispersion-based approach is displayed 
where the quantized conductance steps are observed. In the 
center panel, state-based conduction results are presented. 
Since the Dirac delta-originated peak values depend on the 
broadening parameters, only their relative amplitudes and 
sparsity need to be considered. The energy spectrum un¬ 
derlying the state-based approach is given on the bottom 
panel. All these results corroborate the isosurface analysis 
of Fig. 3. Both the dispersion- and state-based approaches 
pinpoint to the same energy level as the BTE, which is the 
energy level of CS 7 of Fig. 3. It needs to be mentioned that 
in CNWs we have spotted localized states also above the 
BTE . 48 
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FIG. 4. The calculated ballistic electron conductances and energy 
spectrum of 1.5 nm equal-diameter CNW. Top: The conductance 
through the ID channel of [110] NW via the dispersion-based ap¬ 
proach. Center: The transport along the [110] direction within state- 
based approach. Bottom: The energy spectra belonging to the center 
panel. 


Another important remark is that the BTE for CNW is 
very close to the conduction band minimum (and BTE) of 
an isolated NW of diameter of 1.5nm , 35 which affirms that 
indeed the geometric discontinuity of the junction region 
in CNW is responsible for the additional localized states 
below BTE. These findings point out that the quantum 
interference effects play a prominent role at the crossing 
regions. Any discontinuity over the core material silicon 
not only excites higher order evanscent modes 26,52 but also 
results in reflections depending on the state energy and 
leads to standing waves. 

This situation gains a further dimension in the case of 



FIG. 5. Isosurfaces for charge densities corresponding to the low- 
lying conduction band states for a (110) three-wire crossing with wire 
diameters of 1 nm, 1.5 nm, and 2 nm. 


CNWs with all-different diameters. In Fig. 5, we present 
the isosurfaces for charge densities of a three-wire crossing 
with diameters of 1 nm, 1.5 nm, and 2 nm. We observe that 
as the difference among the diameters of crossing wires be¬ 
come larger, the electron states at the band edge delocalize, 
in particular, along the thickest NW. Here, similar to the 
equal-diameter CNW of Fig. 3, the energy of the first ex¬ 
tended state CS 4 nearly coincides with the LUMO energy 
of an isolated 2 nm Si NW . 35 In Fig. 6 , the isosurfaces are 
given for the case of wire diameters of 1 nm, 1 nm, and 
3 nm. Concerning the LUMO, beyond a diameter contrast, 
the CNW behaves as isolated NW arrays, as if the other 
thinner wires are not present. For the higher-lying states, 
while still extending along the thickest wire, they get more 
sensitive to the thinner branches. These figures show that 
the formation and profile of a localized/extended state de¬ 
pends on the relative diameters of the NWs. 



FIG. 6. Isosurfaces for charge densities corresponding to the low- 
lying conduction band states for a (110) three-wire crossing with wire 
diameters of 1 nm, 1 nm, and 3 nm. Carrier densities extend along 
the widest crossing wire. 


B. Kinked Nanowires 

Next, we examine KNWs, here one can benefit from var¬ 
ious opportunities for the electronic structure engineering 
by tuning the diameter and lengths of NW segments or by 
changing the wire crystallographic orientation. Fig. 7 illus¬ 
trates the schematic view of a U-shaped KNW where the 
blue (red) segments are aligned in the [ 110 ] ([ 001 ]) direc¬ 
tion. On the same figure, we specify the variation of LUMO 
energies as the length of the [ 001 ]-aligned segment is de¬ 
creased. As observed from this figure LUMO is particularly 
localized to these sections of the KNW forming a standing 
wave along the [001] wires. Therefore, in agreement with 
the quantum size effect, the LUMO energy increases as the 
length of [ 001 ] segments decrease. 

Focusing on the KNW-3 sample of Fig. 7, on the top 
panel of Fig. 8 , we first subject the KNW to the dispersion- 
based approach which reproduces the conductance steps, 
i.e., NW-like dispersive states, starting from the BTE about 
0.8 eV above LUMO. On the center panel, we present the 
conductance values calculated by state-based approach. 
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FIG. 7. Top: A representative KNW depicted from two aspects. The 
blue segments are along the [110] direction and the red ones are along 
[001] direction. Bottom: Isosurfaces for charge densities of LUMOs of 
three different shape KNWs. Their relative energies are also included. 


Note that the onset of the major peaks is very close to 
the BTE indicated on the top panel. 

We can classify the states in Fig. 8 according to the 
relative values of their conductance peaks: (SI) totally 
localized-bound states, not contributing to the conduc¬ 
tance, (S 2 ) quasi-bound states, having low conductance 
peaks and energetically closer to BTE than to LUMO, and 
(S3) NW-like extended states which have higher conduc¬ 
tance peaks and lie beyond the BTE. As a matter of fact 
these S3 states are responsible for the full ID dispersion 
and contribute to the conductance steps. Dispersion-based 
approach is not able to capture SI and S 2 states, but identi¬ 
fies the S3 states. The validity of this classification becomes 
more founded by correlating these observations with the 
isosurface charge density profiles as presented in Fig. 9 for 
a selection of conduction states. Here, the first 12 states 
(CS 1 -CS 12 ) belong to Sl-type states, indicating localized 
characteristics. The representative states CS 21 , CS 31 , and 



FIG. 9. Isosurfaces for charge densities corresponding to several con¬ 
duction states (CS) of a KNW. 


CS 43 are S 2 -type also marked on Fig. 8 . As observed in 
Fig. 9, their isosurfaces also display an extended behavior. 
Rest of the states given in Fig. 9 are lying at or beyond the 
BTE, exhibiting NW-like extensions (S3-type). 

IV. CONCLUSIONS 



In conclusion, combining wave function analysis with the 
ballistic transport results we map out the localization be¬ 
havior of silicon CNWs and KNWs. We find that the 
BTE can be up to hundreds of millielectron volts above 
the LUMO, filled in between with a substantial number 
of localized states. Their 3D wave functions unambigiously 
display localization around the NW geometric discontinuity 
regions, such as junctions and bends which excite higher- 
order evanscent modes and incubate probability amplitude 
interferences. In contrast to the states close to LUMO, 
those in the vicinity of BTE possess a proto-conducting 
character which might contribute to transport through an 
assisted process, like phonons. Thus, via a realistic atom¬ 
istic electronic structure analysis, we show that ultraclean 
Si NW architectures such as CNWs and KNWs should man¬ 
ifest a conduction band tail that exemplifies a disorder-free 
localization. 


FIG. 8. The calculated ballistic electron conductances and energy 
spectrum of a 1.5 nm diameter KNW. Top: Conductance steps using 
the dispersion-based approach with the corresponding dispersion be¬ 
ing given in the inset. Center: The transport along the [110] direction 
within the state-based approach, also zoomed in the inset. Bottom: 
The energy spectra belonging to the center panel. 
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APPENDIX A: LCBB FORMALISM 

In the LCBB method, 33,34,48 the j th state of the nanos¬ 
tructure can be expressed as the linear combination of bulk 
Bloch bands 


^•(f) = 


-L V C a , . u a Ar)e' k ' 

yCy n,k,j n,k v ' 

n,k,cr 


ik • f 


where N is the number of bulk primitive unit cells within 
the large supercell of the nanostructure, n is the bulk band 
index, k is the wave vector sampling the BZ of the underly¬ 
ing lattice, and a labels the constituent materials, i.e., core 
or embedding media. The periodic part u°^ -(r) of the bulk 
Bloch function can be expanded by the plane waves as 


N g 


n,k 


(r) = 


-= V B a AG)e iG ' f 
nfeV ; 


(n'k'a'\Vpp\nka) = E B^G 1 )* B^G) 

G.G' 

X E V°"(\k + G-k'-G'\ 2 ) 

<j" ,a 

x $ ) e-i$+G-k'-G'yd°: 

Here and are the Fourier transformations of 

atomic pseudopotentials and weight functions, respectively. 

In the case of dispersion-based approach, 48 the separa¬ 
tion of the k -space into k± and fcy (illustrated in Fig. 2) 
leads to a separated eigenvalue problem of the form of 

E (n'k'Jna'lT + Vpp\nk±k\\<j) 

n,k±,<j 

=EAh) E \nkj l{ a), 


where Nq is the number of reciprocal lattice vectors G 
within an energy cut-off, and Q 0 is the volume of the prim¬ 
itive cell. These bulk expansion coefficients can be 

calculated by diagonalizing the bulk Hamiltonian for each 
k point. In our computational implementation, the bulk 
coefficients are computed by employing the empirical pseu¬ 
dopotential method (EPM). 53 

The total single-particle Hamiltonian of a nanostructure 
with the kinetic energy and the ionic potential parts, where 
the latter describes the atomistic environment within the 
pseudopotential framework, is expressed by 

H = T + V PP = -^-+ E W£(Ri)va(r — Ri - <%), 

cr,Ri ,a 

where m 0 is the free electron mass, the direct lattice vector 
Ri gives the position of the primitive cell and d £ speci¬ 
fies the relative coordinate of a particular atom within the 
primitive cell. The weight function W°{Ri) keeps the in¬ 
formation about the composition of the system by taking 
values 0 or 1 depending on the type of the atom located 
at the position Ri + . v° is the local screened spherical 

atomic pseudopotential of atom a of the material a. 

Hamiltonian matrix elements are evaluated within the 
basis set {| nka)} that is {(r]nka) = -(r)} in the position 

representation. The resulting Hamiltonian is diagonalized 
to yield C a t . coefficients. The corresponding generalized 
eigenvalue problem is 

Y^{n'k'a'\f+V P p\nka) C^^E^C^/n'k 1 o'\nko ), 

n,k,cr n,k,cr 

where the matrix elements are 

(n'k'a'lnka) = 5^, E ^(^*^(<5), 

G 


(n'kw\t\nk„) = s u , ■£ y 01 , 


which is solved at definite fey values to obtain the energy 

dispersion Ej(k\\). This form does not put any restriction 
on the treatment of matrix elements. One should only need 
to form a planar fc-grid at a specific fey point. 48 

APPENDIX B: OXIDE PASSIVATION AND ATOMIC 
PSEUDOPOTENTIALS 

An important aspect of the electronic structure calcu¬ 
lations is the surface passivation. The surface of an un¬ 
passivated nanostructure consists of dangling bonds which 
introduces surface states lying in the otherwise forbidden 
band gap. Surface passivation removes these surface states. 
In our implementation, the surface passivation is provided 
by embedding the NW structures into a host matrix of 
SiC> 2 . However, the EPM calculations involving oxygen is 
rather nontrivial and, moreover, the different lattice struc¬ 
ture of Si02, with reference to Si, causes strain effects. 51 
To overcome these obstacles, we introduce an artificial 
monoatomic wide band gap material that computationally 
substitutes silica. After constructing the Si NW-based core, 
all the remaining crystal points within the computational 
supercell are filled with the artificial oxide atoms. This ar¬ 
tificial embedding material has the same dielectric constant 
of SiC >2 and band edge lineup, with respect to Si, however 
with the diamond structure to prevent complications asso¬ 
ciated with strain 51 . 

For the local empirical pseudopotentials of Si and Si02, 
we use the functional form introduced by Freidel et al. 54 
These authors suggest an analytical expression to produce 
the pseudopotential form factor at a given general wave 
number q 


Vpp(q) = 


«i(g 2 - a- 2 ) 

e a 3 (q 2 -a 4 ) _|_ \ 



1 

2 


( 2 ) 


We employ their parameters for Si (listed in Table I). 54 
On the other hand, we have generated pseudopotential 
parameters 48 for the artificial Si02 (given in in Table I) 
which reproduces the experimental band alignments of bulk 
Si/Si02 interface as 4.4 eV and 3.4 eV for the valence 55 
and conduction 56 band edges, respectively. The units are 
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TABLE I. Parameters of the pseudopotential form factors (see 
Eq. (2)) of Si and the wide band-gap matrix which represents SiC> 2 . 
Units are given in the text. 



a± 

CL2 

a 3 

<24 

a 5 

<26 

Si 

106.0686 

2.2278 

0.606 

-1.972 

5.0 

0.3 

Matrix-Si 

69.625 

2.614 

0.0786 -19.1433 

5.99 

0.335 


arranged such that the pseudopotential form factors come 
out in Rydbergs and the wave number q in Eq. (2) should 
be in the atomic units (1/Bohr radius). For the plane-wave 
cut-off energies of EPM, we use 14 Ry for both Si and arti¬ 
ficial Si 02 - 

Though the surface chemistry of oxide-passivation re¬ 
quires more elaborate first principles calculations, consid¬ 
ering today’s computation limits, we believe that our ap¬ 
proach offers a viable alternative for the calculations of 
oxide-embedded large scale Si nanostructures. Despite 
the missing surface relaxation and strain effects, the re¬ 
liability and competence of our computation method has 
gained confidence in the context of embedded Si and Ge 
nanocrystals, 51 and in the energetics of single Si NWs and 
NW networks. 35 

The LCBB basis set is constructed by only employing the 
lowest two bulk conduction bands of the core Si. We tested 
and observed that the further increase in the number of 
core bulk bands or the inclusion of the bulk bands of oxide- 
matrix do not alter the energetics and localization behavior 
of the states within the energy window very close to band 
edge (~0.5 to 1 eV above the conduction band minimum). 
Further discussion of LCBB method can be found in the 
original papers 33,34 and in our previous works. 35,51 
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